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SDBJECTj biweekly report, NOVEMBER 2, 1953 

Toj Jay W» Forrester 

Fromj Scientific and Engineering Computation Group 

!• MATHBHATIGS. CQDIKG AHD APPLICATIONS 
1«1 Introduction 

During the period covered tsy this report 213 coded programs were run on 
the tirae allocated to the Scientific and Engineering Computation(S&EC) Group* 
These programs represent part of the work that has been carried on in '20 ©f the 
problems that have been accepted by the S&EC Group* Progress on each of these 
problems is given below in terms of programming hours j minutes of computer time 
and progress reports as subiiitted hy the programmers in question* 

A new format has been introduced £or presenting these progress reports* It 
is hoped that a more conplete description will be provided thereby for each prob- 
lem* The significant changes in the format include the following: 

(a) The problem ntambers may now be followed "by an upper case letter A, B, G 
or D* The absence of such a letter indicates that the problem originated within 
and is being solved by the S&EC Group* The letters themselves denote whether a 
problem is being carried out for academic credit (e«g,, thesis work) and whether 
the problem is sponsored, that is» whether some organization (outside of M*I*T.) 
is contributing to the financial support of the problem* The letter A indicates 
that a problem is neither for academic credit nor sponsored, B indicates that it 
is just for academic credit, C indicates that it is sponsored but not for academic 
creditj and D indicates that it is both for academic credit and sponsored* 

(b) The titles of the problems have been expanded into abstracts to enable 
the reader to determine the physical and mathematical contents of each problem^ 
These abstracts will be repeated in each biweekly report* Programmers will con- 
tinue to write their biweekly reports describing the progress during the period 
covered (with the exception of their initial reports which should describe the 
nature of the problem, the numerical procedure to be used, and the resuLts de~ 
sired) * 

(c) On the line following the abstract is the name of the person (or persons) 
supervising the problem or directing the project supporting the problem* The 
laboratory or institute department is also indicated* In future reports it is 
planned to include on this line the specific prefect numbers* 
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(d) On the folloTrjing line are listed the names of the persons actually 
carrying out the work on WWI« The name of the author of the biweekly report 
sutraitted for the probljem has been underlined. In problems (e«g» #100 and #140) 
where more than one person has contributed a biweekly report, the author's name 
is inserted at the end of his report. The programmer's laboratory or department 
affiliation will be indicated only when it differs f rem that of his supervisor. 
The following abbreviations iiave been introduced to describe the programmer's 
staff positions 

D«I«G« ; Division of Industrial Cooperation 

CMMC. J Committee on Machine Methods of Computation 

Res« Assist, 5 Research Assistant 

The time indicated for each programmer has been reported by him as the time he 
has spent during the biweekly period on the given problem, 

(e) On the line beginning with DCL (Digital Computer Laboratory) will be 
found the names of staff members of the S&EC Group who have spent the indicated 
amount of time on the problem. Also on this line will be found the amount of 
"WWI computer time constimed by the problem during the period in question* 

1^2 Programs and Ccmputer Operation 

100, Cdm-prehensive System of Service Routines « developed Isy the S&EC Group at the 
Digital Computer Laboratory for the input conversion of suitably prepared 
punched paper tapes, When so requested p these routines automatically 
provide a program with suitable programmed arithmetic, cycle counting, 
and output facilities, 

sIXJL Staff 9 CcMbeliCj 1 hour;Denman, 56 hours; Demurjian, 9 •S hours; 

Frankovichp 20 hours? Helwigj 25 hours 1 Hazel » 10,5 hours; Kopley, ? 

hours; Porter* 29 hours; liMl^ I98 minutes 

The comprehensive system has been modified so that the basic conversion 
program will be recorded on unit preceding the first pass of the CS, The 
input program has been modified to provide two entry points for basic conver- 
sion (one using the photoelectric tape reader and the other using the mechanical 
reader). Using either of these entry points will result in the basic conversion 
program being read from unit into storage. These modifications have been 
successfully tested and will be incorporated into the various programs. 

An error in a modification to the read-in program has been discovered. 
This caused an initial caO in a 356 tape to be treated as an illegal character. 
This mistake will be corrected as part of the preceding modification of the 
read-in program. 

Consideration is being given to the rewriting of the programmed arith- 
metic section of the comprehensive system, 

Helwig 
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The comprehensive system (CS) is being modified so that the displ^ seope 
will no longer be used during conversion* The title of a tape being converted 
will be recorded on the delayed output equipment along with the eonverrted pro- 
gram. If no unassigned flads (floating addresses) are used in the program, then 
the table of flads will also be recorded on this equipiment» However, if aiQr 
unassigned flads are requested,, a table of these will also be printed* 

An error in the generalized decimal number conversion has been corrected* 

All the changes described above wi|Ll take effect the next time GS is re- 
corded on magnetic tape unit 0« 

Frankovich 

The flad table is now being printed out on the delayed printer instead of 
being displayed on the scope« This means that a flad table will be available 
at the same time the 556 tape is made which will save the programmer unnecessary 
delay, Also^ the table has been changed to give the absolute addresses in the 
decimal instead of the octal number system^ 

The program for the generalized decimal number post~mort®n has been moved 
so that it is now located at the end of the second bank of high-speed storage* 
The prograircaer may now have ang?" G»I5?? n"umbers in the range 40(32) through 3203 
(1667) printed out. 

Hazel 

The possibility of using the character generator for output and post-mortems 
is being explored* Results thus fa,T have been fairly satisfactoryo The system, 
as it exists is not particularly well suited for using both the character gener- 
ator and the point=wise display without recalibrationo 

The present set-upj, though far from idealg may have some limited uses© 
However 9 one can readily envision a future development of the system which would 
permit switching from one system to^ the other merely by executing the appropriate 
output instructiono 

Kopley 

Eight chapters of CS Manual I have been written*. The first four chapters 
have been revised^ eorrectedj, and are now available for limited distribution » 
Chapters 5 through 8 are in the process of being corrected and revised, 

CS Manual I is designed to provide an introduction both in programming and 
in the logic of the Comprehensive System* 

Porter 

101 C» Optical Properties of Thin Metal Films on transparent backings are deter~ 
mined and printed out automatically by this program 5 the input data 
consist of the observed reflection and transmission c ©efficients , the 
index of the backing ^ the wavelength and the sample thickness o The 
program calculates by means of an iterative procedure and prints out 
the index of refraction and the absorption coefficient of the filnig 
the rate of variation of these constants wi-tti reflection and trans^ 
missions and the fume's conductivity and dielectric constant «~ 
tfor Professor L<, Harris 9 Chemistry Department 
sby Pro Ao Lj^^iioeb (DIG) „ 15 hours? Jo Richmond(DIC)9 20 hours 
sDGL Staff, Demaano 1^5 howrsi l^MIg I30 minutes 
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During the past biweekly period the eraphasis has been placed on r\anning 
the very thinnest films (thicknesg less than 100 A,®), which give alarms « 
Heavier films h^ve^ with one exceptions given reasonable results without diffic- 
•altyo As. the thickness of the films decreases the optical constants become 
indeterminate* The alarms have been analyzed and traced to this indeterminacy* 
Thus the scope post-mortems provided very useful results for a mathematical analysis 
of the accuracy afforded ty the experimental measurements made on the very thin- 
nest films. 

On the heavier films data published by Fo Goos were used^ and the results 
frcra WWI compared with Goos' own calculationo Goos only made the calculation 
for part of his results ^ and so we filled in the gaps. On the samples for which 
Goos did calculate constants he only calculated one set of constants » while we 
were able to utilize front and back reflection to get separate sets. Goos' 
experimental results give quite a discrepancy for the optical constants obtained 
from, front and from back reflections this indicates that his backings were not 
properly cleaned. 

106 0)4 -* MIT Seis mic Pro^^ect is concerned with the development of methods for 

locating deep reflections from underground strata in seismic pros- 
pecting* The basic method is one of prediction by means of an opti« 
mum linear operator. 

sfor Professor PoM„Hurley<, Geology and Geophysics? Profess or G.Wadsworth, 

Mathematics Department 

?by E.A .Robinson ( Res « AssoCo)oH.Brisc0e9l8 hours? S. Simpson^ 15 lioursi 

W. Walsh p 15 hours 

sDGL^ WIp 266 minutes 

In production runs this week 13 error curves and 12 spectra were computed, 
and 19 new matrices were set up from the cross products of readings from seismic 
x>races * 

A program which uses the delayed printer to produce a continuous plot of 
a long list of positive numbers has been written and tested. The ordinate is 
represented by the number of spaces from the left hand margins and the abscissa 
changes by a fixed interval,^ determined by a carriage return ^ between each point. 
The values of the numbK:^ are also printed outj, as decimal fractions,, along the 
right side of the page, 

A program to check the change in the average of the readings of a seismic 
trace going down the trace was used on all traces that have been taped. This 
gives a check on the reading of the traces since the average should be constant. 

Three attempts to solve a matrix using a program now in the subroutine 
library were unsuccessful^ twice due to tape errors and the third time for an 
undetermined reason. 

108 ©. An Interp retive Program is being developed that will accept algebraic 
equations p differential equations ^ etc* expressed on Flexowriter 
punched paper tape in ordinary mathematical notation (within certain 
lijftits imposed by the Flexowriter) as input and automatically provide 
the desired solution « 
?for Dr. J.HoLaning Jr.^ Instrumentation Laboratory » 
sby JoHoLaning Jr.(DIC),^ hours? N.. Zieri^r: (DIG), 40 hours 
s:DCL Staff 9 Hazel p .5 hours fWWIo 19 minutes 
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The program now carries out the automatic computation of special functions 
satisfactorily a J«H,Laning9 Ji'^o has written routines to compute logarithms, 
exponentials, the hyperbolic functions and a routine for the automatic solution 
of differential equations* These routines will be incorporated into the main 
program in the near future » 

112 C3J* Lawley^s Method of Factor Analysis is applied to a correlation matrix 

obtained frcxn psychological tests and grades given at the Naval Academy^ 
This requires the solution of a modified eigenvector equation for this; 
matrix and is accomplished by an iteration of Hotelling's method for 
solving an eigenvector problem© 

sfor Dr* F« M«, Lord^ Educational Testing Service » Princeton, N«J, 

sDCL Staffs Denman,2 hours i WWI^ 6 minutes 

After correcting the scale factoring of certain numbers, the matrix 
problem was run again with an assumed rani? of 8 for the correlation matrix* Ho 
alarms occurred and these results are now being examined at the Educational Test- 
ing Service* 

119 B* Spherical Wave Fropakation produced ty the sudden release of a spherical 
distribution of compressed air in the atmosphere is being studied by 
numerical means* This involves replacing a set of non-linear hyper- 
bolic partial differential equations in 2 independent and 2 dependent 
variables by a set of difference equations written along characteristics* 
An iterative procedure is used to solve these equations* 

sfor Professor CoC*Lin9 Mathematics Department i, 

^■^ A* Ralston ., CMMC,, 5 hours 

sDGL, WWI, 52 minutes 

Work on the problem as originally stated is now complete* Physical checks 
applied to the results indicate a high-degree of accuracy* For the initial 
density distribution used, the expected shock wave did not materialize* On the 
contrary, after one outward wave from the center and one inward wave to the 
center the state of the system was near equilibrium and there was every indicat- 
ion that wlaat disturbances were left would quickly die out* 

In the future it is planned to try different initial density distributions 
in an attempt to read a shock* 

126 G» A Data Reduction Program for use in the Servomechanisms Laboratory is 
being developed in separate stages to be combined at a later dat6* 
The first stage is concerned with djbvising a program to fit polynomials 
to arbitrary empirical functions using a least squared error criterion* 
The procedure makes use of Legendre polynomials and matrix multipli- 
cation* 
si;or J *E, Ward J Servomechanisms Laboratory;, DIG Ho* 7138 
%hy B* Ross (DIG) 6Q hours? D«Hamilton(DIC), 66 hours i R.Turyn(DIC)80 hours 
sDGL Staff, Frankovich, *5 houri WWI,, 202 minutes 

During this period substantial revisions of the Mistake Diagnosis Routine 
(MDR) were made to correct two logical errors and reduce the number of restrict^ 
ions on which orders can be replaced by break points* Test runs have checked 
all but one section of the program and have disclosed some pitfalls in the 
application of the MDR* Since the MDR now will satis ly the needs of most users » 
problem 126 includedg no further developmental work on the MDR is contemplated* 
A report describing fully the application of the MDR is almost finished and will 



be published as Serromechanisms Laboratory Report 713S~SR'=-1 mtiiout delay« Copies 
■will be available from the laSiittemore library® Persons -Mho desire copies from 
the initial distribution may call extension 2liil and leave their name and room 
number* 

Mi. program for Lagrange Interpolation using the n-point^ equal spacing for- 
mula has been written^ The program is -written in terms of preset parameters which 
specify the number of points (n)^ the number of distinct functions which are to 
be interpolated at the value of their common independent variable which is in 
the ySSA upon entry to the routine ^ and the pattern in which the given functional 
values are stored© fhe program may be considered to interpolate values in an 
m-dimensional vector function^ Initial tests have not been successful and reasons 
are being soiight* 

4 post-mortem routine to printc, as generalized decimal numberS|, arbitrary 
blocks of (21^^ 6) numbers stored anyidiere in the coiiputer ©r on the drum has been 
written and is ready for conversion^ The locations of the blocks are specified 
by an easily prepared Flexo tape to be read by the mechanical tape reader either 
by initial and final address (oc°tal or decimal) or by initial address (octal or 
decimal) and block length (decimal only)*, Drma groups are specified decimallye. 
Ho Benman^s (2i|.5)6) print routine is being used for oi 



^^^ ^^ §!^2^^^^^ SSF^^R.^^^nr^^^'^XS^^^^^^^ ^°^^=^JS£iS;B2 ^^^ being revised^ 

^nd~^Eested« "The resultii^^ set of rouCineF^an be used when desired to 

obtain the necessary data^ The subroutines involve routine numerical 

and logical operation,?* 
gfor Professor ¥J'L Pease ^ Servomechanisms Laboratory 
sby J»HoRun3^ri (S.si*Res. Assistol^JS hoursj 
§DGL .Staffs E^ankovich^ 1 hour^- Wl^ 90 minutes: 

Two more library subroutines were successfully tested* This conpletes the 
testing of the basic set of subroutines <> Three subroutines for finding points 
on the series l6 symmetrical wing section^ camber linei, and asymmetrical wing 
section were written and are being tested©. 

& rather baffling error in one of the routines successfully tested was: 
traced to a change in the s^lcx®- order idiich had not been noted© 

A tape for the supersonic nomzlB was obtained© This tape has 2530 milliaag 
machine blocks* The run which produced it took 1? mnutes of -sdiich approximately 
five minutes were required for output using delayed printing and punching* The 
punching of the tape requires about forty minutes ©. This tape is for a k5° segment 
of the nozzle© Since the nozzle is formed by rotating a given cross section 
through 360^0 the entire nozzle can be obtained from the kS^ segment by inter** 
changing millii^ machine axes a Time spent in aetual cutting is expected to b© 
approximately 26 hourSe^ 

136 Matrix Equations© Yarious metoods have been studied for the solution of 
a. set of lanear algebraic equations^^ A variation of the Hestenes— 
Stief el conjugate gradient method has been prograssaed and tested for 
insertion In the SSeEC Library of Subroutines 
§DGL Staff I, ilrden^ 2: hours| Helwig^ 10 hourS| Wl<, 4- minutes; 

Tape 262iimll will be included in the WH library of sub3?outines«^ This 
program will solve a set of linear equations which are stored on the drum in 
the form of an augmented matrix in column order ^ ie©*^ the first column of 
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coefficients and so on for n colmns of coefficients ("where n is the maiaber of 
equations in the system) and finally the coltjuin of the non-homogeneous parts of 
the equations. Two preset parameters must be specified! 

zml = 2n 

zm2 *= initial drum address of stored data 

The program itself occmpies li^5 registers , and 13 + ^saal registers of 
temporary storage are required, "^he original matrix is unaltered Toy the program • 
The program operates on an initial guess stored in zml registers beginning at 
13t + zml and calculates an approximate solution, which is stored in zml registers 
beginning at 13t + zml. %e error resulting from the substitution of this solution 
in the original equations is also calculated and is stored in zml registers begin- 
ning at 13t +2zml. This approximate solution can very often be improved by simply 
re-entering the routine without altering the contents of the zml registers begin^ 
ning at 13t + zml» 

The program must be entered and is left in the interpretive mode. 

137 B, Investigation of Atmospheric Turbialenee as a noise input to airborne 

control systems is being studied as a stationary random process. 
Hence the methods of generalized harmonic analysis may be used to 
describe the turbulence in terms of its power spectral density. 

sfor Professor R.C.Seamans, Aeronautical Engineering, Instriimentation Lab. 

sbv Rf A. Summers . 10 hours; N. Zierler, 5 hours? G. Block, 4 hours; 

:DCL, Wl, 160 minutes 

The lOO-aOGtr autocorrelation functions for runs I, II, III, and 17 were 
jsompleted using data tapes clearly marked to indicate the fifth block of data. 
Agreement with the 0-100 r functions was satisfactory for runs II, III, and IV, 
but not run I. The reason for the discrepancy in run I is not known at the 
present time. 

I^ the Fourier transform program of problem 107, an error was foui^d in 
T-2235p9 and corrected in T-3263pO» but the apparently spurious negative' bias 
in the power spectra is still observed. In order to detennine whether this 
bias really occurs in the data or is due to a program error, N. Zierler has 
utilized the generalized program of problem 108 to compute the Fourier transform. 
This is a program cf considerably higher accuracy, but, due to the increased 
machine time required* it will probably be used only once or twice as a check on 
the more rapid transform program of problem 107» The Fourier transform has not 
yet been successfully computed ty the method of problem 108. 

Completion of problem I37 now depends on satisfactory operation of the 
Fourier Transform program. 

138 B. Spheroidal l^^re Functions are solutions of the scalar Helmholtz equation 

separated in spheroidal coordinates. A program has been developed for 
tabulating the solutions of the first kind in terms of the coefficient's 
of their exptosion in associated Legendre functions. These coefficients 
are determined ty applying an iterative procedure to the difference 
equation approximation of the corresponding separated ordinary differ- 
ential equation. 

sfor Professor P.M.Morse, Physics Department 

sby F. J .Gorbat6 ,( Res. Assist. CMMG),> 100 hourssJ.C.Little(Res.Assist.GMHG) 

jDCL Staff, Combelic, 3*5 hours §1M, 83 minutes 
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it coK5)lete prograin for tjiis problem has iiow been tested successfully* l&e 
program has already been used to investigate the extreme values of the parameters 
and the 'results indicate a possible minor modification in the normalisiation 
conditions© This theoretical consideration -will be settled shortly and then full 
production is anticipated* 

!2he following significant difficulties were encountereds 

(1) ii^en a programmer desires to use more than one buffer and therefore ^lses 
the auxiliary buffer subroutine^ the P& section of CS is altered so that it 
refers to this subroutine* Consequently if this buffer subroutine is read over 
by another routine (e*g»3, by exchanging it teit^jorarily mth a routine that has 
been stored on the drum)^ the PA may not operate satisfactorily* To avoid this 
difficulty the programmer should regard the auxiliary buffer subroutine as an 
integral part of the Pit* 

(2) During conversion^ a *»ca 0*«' is punched out every I4.8 feet of pxmched tape and, 
if by chance this falls after an wfep block^^ the *»tea 0® is erroneously interpreted 
as an illegal first character of the next block* . The correction needed to avoid 
this difficulty idll be incorporated into the read-in program (cf ♦ problaa #1(30)*. 

lliO* Summer Session System consists of a conversion program, on interpretive 
routine, and mistake diagnostic routines stored in IMl* ^ special 
mnemonic instruction code has been developed for "use "Hit^ this system 
thus simulating a coii^juter "Bsith characteristics quite different from 
those of Wl* This Summer Session (SS) cosputer was developed for the 
use of students participating in the MET 1953 summer session course on 
"Digital Cbi!5)uters and Their Applications^^'* The IS coBfJuter is being 
used in the E*E* Department course 6^535 and is available 1^ programmers 
with suitable problems* 
sDCTi Staff, Siegel j) SO hoursj Combelic, 50 hours| Denman, 5 lioursi Frankovlch, 
IS hours| Helwig, 25 hours 1; Hoy, 59«5 hours | ^JSl, S^ minutes 

The siaamer session (SS) coi^uter was used by the stiadents in Course 6*535 
for a five-feour computation period* The operation of the computer was satisfact- 
ory in every respect* Copies were Retained of all postHsiortems produced during 
the coB^utation for a study of the types of mistg^es made by the novice prograi»» 
mers* 

The program for recording the *«SS coii5)uter»* on a magnetic tape unit and for 
reading it directly from magnetic tape to the dnjm has been wcitten and tested* 
It is e3q)ected that the BS coB5)uter "wiH soon be permanently recorded on Unit ^ 

A subroutine is being written for autoraiatically recording the title ©f the 
program on those output media iftiich the programmer has selected before the per- 
formance of the program* This will permit positive identification of typed 
results and oscilloscope displays obtained using the computer* 

The exfei3:p4tion of minor and obscure mistakes in the coir^juter subroutines 
continues* 

iM^ MBS Subroutine Study has been undertaken for the final testing of sub- 

routines selected for incorporation into the Library of Subroutines* 
Although very little effort is going into the specific development of 
subroutines, programs that have been written for other S&IC problems 
and seem to be of general use id.ll be suitably modified for "the S&0J 
Library* 
^K3L Staff, Arden^ 23^ hoursf Wl, 33 minutes 
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A subroutine for solving systems of linear equations by a method developed 
and described under problem #136 has been written and successfully tested* The 
subroutine will be incorporate?! in the WWI library • 

A subroutine adding a format to subroutine tape 2756i5 has not yet run 
successfoolly* 

142 D, A Study of Shock Waves has been undertaken in two dimensional solids 

subjected to impulsive loads. The analysis approximates the solid 
by a two dimensional grid with concentrated masses at nodal points* 
The response of this system is computed from a finite difference 
approximation to the differential equations of motion of this system* 
:for Professor C.H.Norris, Department of Civil and Sanitary Engineering 
sby R* Bart(DIC),60 hours? S> Sydney (Res* Assist* CMMC), 60 hours 
sDCL Staff, Kopley, 1 hour; WWI, 1^36 minutes 

Final results for the analysis of a two-dimensional solid under the action 
of a parabolically distrtbi;ited dynamic load have been obtained* The cooperation 
and assistance given by members of the staff to the programmers during the course 
of this problem is greatly appreciated* A copy of the report on this project 
will be sutanitted to the Digital Computer laboratory by S* Sydney and R* Bart 
of the MIT Department of Civil and Sanitary Engineering* 

1^3 D* The Vibrational Frequency Spectrum of a Copper Crystal is to be determined 
by solving a 3 x 3 secular determinant, each term of which consists of 
a finite Fourier Series of 12 terms* This equation must be solved for 
Z5»^5 different values of the wave-propagation vector* 
sfor Professors B.E.Warren and J ^C .Slater, Physics Department 
sby E.H.Jacobifeen ( Res > Assists )> 8 hours s FoJ.Corbat6(Res. Assist.CMHC) , 
sDCL Staff ,'^'Gombelic, *5 hour; WWI, 16 minutes 

Most of this period has been spent on evaluating a double-intfeeral (for 
27 values of a parameter) which is concerned with the "isecond order diffuse 
scattering" of X-rays. Both integrations involve some thirty values of the 
integrand. The integration uses Simpson' s rule and an inverse tangent routine 
devised by Donn Gombelic. The program is proceeding satisfactorily* 

1^ C* Self -Consistent Molecular Orbitals are the optimum choices of linear 

combinations of atomic orbitals determined through a process described 
as a self -consistent field approximation. The nxamerical procedure involves 
' matrix- vector muliplications , vector additions, and matrix diagonal- 

ization* 
sfor Professor J. C .Slater, Pl^sics Department 

s by Dr. A. Meekler( DIG. Solid State and Molecular Theory Group) ,5 hours 
sDCL Staff, Arden, 13 hours; WWI, 28 minutes 

The program has been held up because of the transfer of a zero to a register 
which is to hold a divisor* Various manual modifications are being used which 
will partition the program in an attempt to learn the time and place of this 
unwanted transfer* 
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1^7 C. Energy Bands in Crystals are being studied tjr finding solutions of the 
corresponding second order linear differential equation satisfying 
boundary conditions at the origin. The solutions are found approx- 
imately by using the Gauss-Jackson formula for forward integration. 
The solutions and their first derivatives are to be combined in a 
sum, the weighting factors being functions of an independent para- 
meter. 

ifor Professor J. C .Slater, Physics Department 

:by Dr. D.J.Kowarth (DlG. Solid State and Molecular Theory Group) , 24 hours 

jDGL Staff, Arden, 9 hours; WWI, 83 minutes 

The routine to calculate the logarithmic derivative of the wave function 
and to repeat the process for a succession of all the parameters has been tested 
and made error free. The final part of the program, to compute a weighted s\am 
of these functions, is being tested. 



1.3 Operating Statistics 

Computer Time 

The following indicates the distribution of WWI time allocated to the 
S&EC Group. 



Programs 


56 hours, 


22 minutes 


Conversion 


10 hours. 


15 minutes 


Hagnetic-Dr\2m Test 




21 minutes 


Scope Calibration 


1 hour , 


21 minutes 


Magnetic Tape Check 




03 minutes 


Total Time Used 


68 hours, 


22 minutes 


Total Time Assigned 


71 hours, 


14 minutes 


Usable time, percentage 


95.95^ 




Number of rrograms 


213 
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2. CGMPUTER EMGIMEERING 

2.1 WWI System Operation 

2#11 Core Memory (N. L» Daggett) 

Now that the Core Memory has been moved to its final position, and a new 
air-conditioning duct installed, the work of putting in permanent video and 
power wiring can begin. The Core Memory Control has already been mounted in rack 
EG with complete marginal-checking coverage available. Transfer of the control 
was 'accomplished very smoothly thanks to the fine job of preparation done by the 
various systems technicians involved. 

(L. L.Holmes, A. J. Roberts) 

The installation of permanent wireways, video cabling, and power wiring in 
the Gore Memory region will be done on regular installation days and will prob- 
ably be completed by November 9» 

2.12 Drum Parity^heck System (L.Loliolmes, A.J.Roberts) 

At the conclusion of the installation day of November 2, the necessary changes 
to the present system to include a parity check of the magnetic dnams will be coni« 
pleted. Room 156 changes might not be done at that time. 

The former Read-In Interlock panel has been converted to the lOR Aux., Mod 
II, and will be located in rack G3. The Parity Auxiliary panel will be replaced 
by a plug"in~unit mounting panel. 

2.13 WWI Service File (D.A<^orrison) 

Another extension to the WWI Service File has been set up. Room 222 has be^n 
supplied with a file catinet to hold drawings pertinent to the equipment located 
in Room 222. Libby Leighton will supervise the filing of the drawing^. Requests 
for drawings in Room 222— «19WI Service File— should be siade to Julie Dickie. 

2.14 Test Programs (S« E. Desjardins) 

The last biweekly period was spent in developing a five- order read-in program 
for use in the five flip-flop registers. This program blocks in a leader tape 
placed in front of the test tape to be read in* The leader tape contains a read- 
in program which then reads in the tape in question. If the regular re«d-in pro- 
gram is used, the leader iape is ignored, and the test program is read in. 

Work was also done on an "inchworm" program combining the inchworm facility 
with a Core Memory test program. More work needs to be done on this. 

2.15 Magnetic Drums (H. L. Ziegler) 

After several days of concentrated effort by Phil White, all relays on one 
Group Selector chassis were induced to operate and release within the required 
time liBiits. This chassis, originally removed from service because of faulty 
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operation » performed quite satisfactorily when returned to service in the 
aiixiliary drum« A second chassis is now being tested and adjusted in the 
Test Rack» 

Because of the maintenance difficulties, poor reliability, and the not-too- 
great life expectancy (conputer-wise) of these rejLays, various alternate methods 
of group selection are being considered. At present the straightforward method 
of individual writers for each group seems most practical* 

(K. £• McYicar) 

The replacement auxiliary drum has been installed in the system and has been 
thoroughly tested • The heads have all been readjusted by an ERA representatives 
and we now get uniform readout signals* Preliminary checking indicates that the 
new copper-clad head pieces do not expand significantly with temperature variations, 
and no attempt is being made to temperature-control the new drum. 

Work on a parity system for the auxiliary drum is proceeding satisfactorily* 
The necessary hardware has been ordered or is being made in the shop* The elect- 
ronic assemblies are now being constructed, and the bay wiring is bein)^ started* 
It has been decided to do group selection of the parity digit by means of crystal 
selection of separate writers instead of using the relay-switching system used 
in the information digits* This decision was the result of the previously report- 
ed trouble we have been having with the relays and the estimated relay life basad 
on latest data concerning programming needs for group switching* 

A review of relay switching for group selection is being made in light of the 
e^^eeted relay life and time consumption with present program requirements* 



2*16 Fair child Camera (L*H«Norcott) 

Recent intermittent failures of the Fairchild Camera to index were apparently 
caused by a faulty tube in the automatic camera-control panel* No failures have 
been logged since the tube was replaced* 

2*2 Terttiinal EouiTMent 

2*21 Ferranti PETR (J*P.Stirman) 

Preliminary tests on the tape reader have been satisfactory so that a bread- 
board is now being constructed for testing in ISWI* 

The brake and clutch circuits are able to stop the tape, from maximum speed, 
in less than one millisecond* Tests *0n the reader circuit have shown that the 
circuit is insensitive to parameter variations within ± 25^ of nominal values* 
These Goi^clusions are derived from performance curves which plot tolerance vs. 



supply voltage. The mlia±mm performance level was set at a signal«otitptit level 
of at least +1 volt and a no-signal output of less than -20 volts. 

2«22 Magnetie-Tape System (E.P«Famsworth) 

The reel of acetate-base tape on unit 3B has broken repeatedly under normal 
operating conditions, and the original mylar reel was replaced. Althoxigh some 
improvement in acetate strength is gained by maintaining it in an atmosphere of 
high humidity, all users of the Raytheon Tape Units seem to have the same tape- 
breakage trouble. Jti is continuing efforts to obtain 1200-foot by 1.5-nJil mylar 
sheets frcia Dupont on a priority basis so that they can supply us with addition- 
al spare reels. 

Programmers making "permanent" recordings on unit are being requested to 
make duplicate recordings on a spare reel of tape on unit 1. Since units and 1 
are now interchangeable by means of the newly-installed locking switch, the exis- 
tence of a duplicate repording for unit 1 will eliminate any inconvenience which 
could be caused by failure of the read/record head or tape on unit 0. 

2.3 Drafting 

2.31 New Multilitji Master,s (A. M. Falcione) 

The attention of all secretaries is called to the new form multilith masters 
which are now available in the stock room. This form is made by Multilith and 
is prei-printed with non-reproducible ink having the proper image area for typing 
reports, memos etc. The dotted image is to be used when a report is being written 
for two-sided printing; this image applies to even-numbered pages. A more detail- 
ed memo willte written to all secretaries regarding this subject in the very near 
future. 
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3« LIBRARY ACCESSIONS LIST 

The follovd.ng material has been receiveci in the Library, W2-325* 
Library Files 
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25^5 


Willow Run Rsch.Ctr. 


2548 


Lincoln Lab« 


2553 


Ballistics Rsch.Lab. 


B>268 

B-269 
B-.270 


John Wiley 

McGraw-Hill 

Prentice-Hall 


Laboratory Files 


No. 


Author 


M«.2J^16 
M=»2425 


G« Young 


M-24ilO 
M«2^1 
M«2459 


R, PfaiTf 


Barta, 


The following material 1 


Library 


Files 


loa^ 


Author 



Title 



A Ccxnparison of Real and Simulated Automobile 
Suspension Systems 

An Application of Matrix Theory to the Tracking 
and Smoothing Equations 

The Propagation of Error in Numerical Integrat- 
ions 

Principles of Transistor Circuits 

Electronic Analog Computers 

Partial Differential Equations in Engineering 
Problems 



Title _>-. 

S&EC Biweekly, 9»20«53 

Digital Techniques for Sorting by Areas in a 

Plane 
S&EC Biweekly for October 5» 1953 
Laboratory Personnel for October 1953 
Marginal Checking for Circuit Designs 



xxtxe 



B 270 KoS .Miller 



Partial Differential Equations in Engineering 
Problems 
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^* AmiNISTRATION AM) PERSOMEL 

Mew Non-Staff ( Ro A» Osborne) 

Tilda Finnochio is a new secretary in the Purchasing Department* 

Roseanne Gillette is a clerk who has transferred from the Bursar* s office to 
work in the Print Room* 

Michael SoH,cmita is a new Laboratory Assistant in Group 63^5 • 

Terminated Non^-Staff 

Vincent Guzziere 
Robert Kyle 
Lois Rutland 
Marlene Wise 



